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Motivated by the striking promise of quantum computation, Majorana bound 
states (MBSs) [1] in solid-state systems [2, 3] have attracted wide attention in 
recent years [4—16]. In particular, the wavefunction localization of MBSs is a 
key feature and crucial for their future implementation as qubits [2, 3]. Here, 
we investigate the spatial and electronic characteristics of topological supercon¬ 
ducting chains of iron atoms on the surface of Pb(llO) by combining scanning 
tunneling microscopy (STM) and atomic force microscopy (AFM). We demon¬ 
strate that the Fe chains are mono-atomic, structured in a linear fashion, and 
exhibit zero-bias conductance peaks at their ends which we interprete as sig¬ 
nature for a Majorana bound state. Spatially resolved conductance maps of 
the atomic chains reveal that the MBSs are well localized at the chain ends (< 
25 nm), with two localization lengths as predicted by theory [17]. Our obser¬ 
vation lends strong support to use MBSs in Fe chains as qubits for quantum 
computing devices. 

Majorana fermions are real solntions of the Dirac eqnation and by definition fermionic 
particles which are their own antiparticles [1]. While intensely searched in particle physics as 
nentrinos, Majorana fermions have recently been predicted to occnr as qnasi-particle bonnd 
states in engineered solid-state systems [2], Snch systems not only offer the possibility to 
observe the exotic properties of snch MBSs, bnt also open np an interesting playgronnd for 
topological qnantnm compnting [2, 3]. The fnndamental ingredients to generate MBSs in 
semicondnctor-snpercondnctor heterostrnctnres is to combine a spin textnre with an s-wave 
snpercondnctor (SC) allowing to create a snpercondncting state with effective p-wave pair¬ 
ing, giving birth to a new state of matter-topological snpercondnctivity [2, 3]. MBSs arise 
as zero-energy states lying in the snpercondncting gap and being spatially localized at the 
interfaces. 

To generate spin textnres, theoretical proposals have snggested to employ nanowires and 
magnetic chains with strong spin-orbit interaction [4-9] or with self-tnning RKKY interac¬ 
tions [10-12], Until now, only few experimental works have reported snccessfnl observations 
of a zero-bias condnctance peak (ZBP) via transport measnrements in semicondnctors in¬ 
terpreted as MBS signatnre [13-16], however, withont addressing their spatial localizations 
in detail. Remarkably, Nadj-Perge et al. [9] recently stndied the spatial and spectral reso- 
Intions of the MBS in Fe chains on snpercondncting Pb by STM and reported the striking 
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FIG. 1: Self-Assembled Fe-Chains on Pb(llO). a, Model of the experiment: topological 
superconducting phase arises when mono-atomic iron chains with spin texture are grown on an 
s-wave superconductor [here Pb(llO)]. MBSs are then localized at the end of the chain and 
experimentally observed as ZBP in the conductance, b, Topographic STM image of iron chains self- 
assembled on atomically-clean Pb(llO). c-d, Topographic STM images of two chain terminations 
with a difference of apparent height (ps 10 pm) at their ends, (V) = -10 mV, It = 100 pA). e, 
dl/dV point-spectra obtained at the chain end in c showing a ZBP. The proximity gap equal to 
A = 1.1 meV was measured with a metallic tip. 

observation of a ZBP at the end of the atomic chains, as generically expected for MBSs [3]. 
Experimentally, the induced proximity gap probed in the Fe chain is found to be very small 

meV), while the exchange interaction is in the eV range, implying on theoretical grounds 
a large localization length of the MBS wavefunction, in contrast to the observation [9] have 
triggered interesting discussions on the physical origin of the ZBP [18, 19]. It further implies 
that possible MBS in such chains may easily hybridize into conventional fermions, raising the 
question as to whether MBS in such Fe/Pb hybrid system will exhibit non-Abelian braiding 
statistics [3]. 

In this work, we go a decisive step further and combine STM and AFM measurements 
at low temperatures to characterize the electronic properties and structure, respectively, of 
superconducting Fe chains self-assembled on the surface of Pb(llO), see Fig. la. To engi¬ 
neer such systems, atomically cleaned surfaces of an s-wave superconductor, Pb(llO), were 
prepared by a few sputtering/annealing cycles (see Supp. Info Fig. SI) on top of which Fe 
atoms were deposited on the sample kept at ~ 400 K. In such way, topographic STM images 
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reveal the formation of long chains extending up to 70 nm (Fig. lb). In agreement with 
[9], chains are partially covered by small clusters mainly localized around their centers (see 
Supp. Info Fig. S2). Careful real-space observations reveal that numerous Pb vacancies are 
also created in the vicinity of the chains. We assume that the relative high temperature re¬ 
quired for the self-assembly also promotes the diffusion of Pb atoms and their nucleation at 
the chain sides. During our study, several chains having lengths from 40 to 70 nm (Fig. lb) 
were investigated and showed two distinct topographic signatures by means of STM at their 
ends (Fig. 1 c and d). They correspond to a variation of the apparent STM height of about 
10 pm. Along the chains and independent of the termination types, no clear atomic period¬ 
icity has been observed by STM (see Fig. S3) which could unambiguously be related to the 
structure of the chain. We then systematically acquired conductance point-spectra at the 
two chain ends to identify the presence of a ZBP which we interprete as signatures for MBS. 
Figure le shows a typical single-point conductance spectrum at the end of the chain with 
end-state (Fig. Ic) obtained with metallic tip. The superconducting gaps measured with 
such tip are thus equal to A = 1.1 meV being the superconducting gap of lead. Therefore, 
we particularly focus on those chains which exhibit end-states, as they show a clear signature 
of zero-energy end states in contrast to chains such as shown in Fig. Id (See details in Supp. 
Infos Fig. S4). We note that the lowest temperature of our experiments is 4.7 K which 
limits our spectral resolution to about 1.4 meV. Although the energy gap in our spectra is 
broadened compared to [9] due to the higher measurement temperature, a superconducting 
proximity gap of about = 1.1 meV is observed in the Fe chain, as well as a clear ZBP which 
we interprete as a signature of a MBS that has been predicted to occur in such systems due 
to various mechanisms [7-12]. 

Electronic and Structural Characterization along the Mono-atomic Fe chain. A ZBP lo¬ 
calized at the end of the chains is one of the hallmarks of a MBS [9]. Alternatively, however, 
such a ZBP can arise from magnetic impurities such as Shiba (near) mid-gap states in the 
vicinity of single adatoms [20, 21], molecules [22, 23] or disorder effectss. To further identify 
the origin of the observed zero-energy subgap states and attribute it to a MBS, we compared 
STM and AFM imaging obtained at the atomic scale below and above the superconducting 
transition temperature (Fig. 2). 

Figure 2b shows a close-up STM topography of the investigated part of a long chain. 
The chain exhibits the topographic signature of an electronic state at its ends as well as 
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FIG. 2: Zero-Bias Electronic and Strnctural Characterization of the Fe Chain, a, STM 

topographic image of the Fe chain end {Vt = -11 mV, It = 80 pA, scalebar = 5 nm). The inset 
shows the section characterized in the rest of the figure, b, STM topography at V = 10 mV and 
/ = 10 pA. c, Constant-height zero-bias AFM image of the chain end acquire at 5 K. In contrast 
to STM topography, the AFM image clearly resolves the atomic structure of the Fe chain. Each 
protrusion corresponds to single Fe atoms aligned in a straight fashion along the Pb(llO) rows 
(Fe inter-atomic distance = 0.37 nm). In c site 1 shows the MBS and 2 the presence of slightly 
misaligned atoms along the chain, d, Zero-bias normalized conductance maps oc LDOS acquired at 
5 K. At site 1 where the MBS is localized, a zero-bias conductance peak is clearly observed at 5 K 
and vanishes at 10 AT as a result of the suppression of the superconductivity. The LDOS obtained 
at 10 K shown in the inset is homogeneous along the chain whereas weak oscillations are observed 
at 5 K due to the MBS wavefunction decay, e, Comparison of LDOS profiles extracted along the 
chain at 5 K (black dots) and 10 K (grey dots). The orange curve shows a fit of the LDOS with 
the theoretical MBS wavefunction composed of two localization lengths. 

a weak perturbation of the electronic density along the chain. This last feature might be 
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caused by the presence of defects below or within the chain. Howewer, since STM reflects 
the electronic density between tip and sample, the ‘true’ atomic structure can be masked 
by delocalized electronic states of the system [24, 25] . As we noticed for the present system, 
the STM topographic data can lead to important misinterpretations of the atomic structure 
of the chain. 

To unambiguously address this issue, we employed the AFM imaging technique which 
is rather insensitive to the delocalization of electronic states close to the Fermi level. Fig¬ 
ure 2c shows a const ant-height zero-bias AFM image obtained at 5 K reflecting the true 
atomic structure of the chain. Each protrusion corresponds to a single Fe atom which are all 
aligned and centered between the atomic rows of the Pb(llO). The chains are mono-atomic 
and strictly ordered in a linear fashion. The Fe inter-atomic distance is equal to 0.37 nm 
(see prohle in Supp. Infos. Fig. S4), which is in good agreement with the atomic lattice 
along the Pb rows (about 0.35 nm). This suggests a high commensurability of the chain 
periodicity with the underlying substrate (the mismatch is about 0.6 %) in contrast to the 
STM observation (see Supp. Infos. Fig. S3). 

The site marked by the arrow 1 in Fig. 2c (X = 2 nm) corresponds the MBS location and 
shows a less negative A/ of about -2 Hz compared to the rest of the chain (A/ background 
= -20 Hz, see prohle Fig. S3). This A/ variation is localized at the last two atoms of the 
chain and corresponds to a less attractive force regime between tip and sample. The AFM 
map (Fig. 2c) shows a perfect round-shape halo of 0.8 nm in diameter as a spatial signature 
of this peculiarity. Within the halo, the last atoms are visible and perfectly aligned with the 
rest of the chain thus excluding atomic disorder or corrugation effects. The decrease of the 
attractive forces at such close tip-sample distances implies additional contributions of repul¬ 
sive forces between tip and sample. Tuning fork AFM is known to be sensitive to such forces 
which usually originate from Pauli exclusion principle between the electronic wavefunctions 
of the tip and the surface [24]. Given that two end MBSs can host at most one fermion, it 
is plausible that a similar repulsive Pauli effect is responsible for the observed halo in the 
force maps. In comparison, site 2 shows a less pronounced A/ variation compared to site 1 
which coincides with slightly misaligned Fe atoms with respect to the chain structure. The 
A/ is less negative by ~ -10 Hz compared to the chain. We further point out the absence of 
any halo at site 2 and conclude that the A/ variation is imputed to atomic disorders, most 
likely due to local topographic variations. We note that we also investigated chains without 
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end-states and no such force contributions were observed (see Supp. Infos Fig. S5). 

To further support our observation of MBS, we compared normalized conductance maps, 

1. e. {dI/dV)/{I/V) oc local density of state (LDOS) between tip and sample at the Fermi 
level, obtained at the same locations and at about 5 K and 10 K, respectively. Fig. 2d shows 
the LDOS{x, y) map at 5 K and reveals a clear ZBP, ascribed to the MBS, localized at the 
end of the chain as determined by the AFM data. A weak modulation of the LDOS is also 
observed along the chain which we attribute to the decay of the MBS wavefunction. At site 

2, the LDOS is almost zero and might be due to weak magnetic perturbations induced by 
lattice disorder [20, 22, 23]. In order to suppress the superconducting state of the system and 
thus to enforce the disappearance of the MBS, we measured the same chain above the critical 
temperature of lead (T > Tc = 7.2 K). The corresponding normalized LDOS{x, y) map (see 
inset Fig. 2d) shows a homogenous LDOS along the chain which drastically contrasts to the 
one at 5 K. By comparing the LDOS prohles taken along the chain at 5 and 10 K (black and 
grey dots in Fig. 2e), the ZBP has completely disappeared at 10 K (site 1) as well as the 
oscillations along the chain due to the suppression of topological superconductivity. Since 
no external magnetic held is applied [9], we clearly address the interplay between supercon¬ 
ductivity and the ZBP observed in our data which provides a strong evidence of the MBS 
presence in this system. 

The MBS Localization Lengths. To further exclude other explanations of the ZBP [26], 
we next investigate another strong hngerprint of a MBS, namely its wavefunction and asso¬ 
ciated localization lengths [17]. For this we assume that the Fe chain with a proximity gap is 
driven into the topological phase by a spin texture which gives rise to a helical held (see Fig. 
la). Although Nadj-Perge et al. [9] measured the chain magnetization and concluded to a 
ferromagnetic behavior, the precise value of the magnetization is however not known and 
could be substantially away from full ferromagnetic order. In turn, this does not exclude 
helical order since the helix can be around the magnetization axis with some small angle. 
Moreover, we wish to point out that, recently, such a helical state has been experimentally 
reported in a similar ID-system [27]. Such helical helds can either follow from spin orbit in¬ 
teraction combined with Zeeman helds [5, 6, 8] or exchange helds [7, 9] or from a ‘self-tuning’ 
RKKY interaction between the Fe spins [10-12]. These two mechanisms are related by a 
spin-dependent gauge transformation and thus mathematically equivalent [28]. We further 
assume that the pitch of the helix is much larger than the lattice constant of the Fe chain 
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FIG. 3: Majorana Localization Lengths, a, Constant-height zero-bias AFM image of the 
defect free chain, b, Zero-bias normalized LDOS of the chain hosting a MBS. c, dI/dV{x) profile 
(red curve) taken along the chain revealing the localization lengths of the zero-bias conductance 
peak (x = 0). The blue curve corresponds to the probability density IV’P of a Majorana Bound 
state with two localization lengths ~ 22 nm and ^2 ~ 0.72 nm. They correspond to 59 and 2 
atomic sites respectively. The orange curve approximate the wave-function by considering the 

effect of the tip radius d = 0.17 nm. 


and that the helical held and the spin orbit interaction is snfhciently strong. In this case, the 
MBS wavefnnction is a snperposition of contribntions coming from different extrema in the 
spectrnm, which resnlts in two different localization lengths for a single MBS [17]. Assnming 
a semi-inhnite chain with one MBS at each end (withont overlap), these localization lengths 
are determined by the corresponding gaps and given by 

hvp 


A ’ A„,-A' 

The MBS probability density is then dehned as 

\^(x)\^ = l{e-2-/«i + e-2x/6 

— cos{2kpx)}, 
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where N is the normalization constant, and ^2 are the localization lengths, A and 
are the proximity gap and the helical held gap of the chain (which depends on the exchange 
coupling constant), respectively, and kp is the Fermi wavevector of the chain. The topolog¬ 
ical phase of the chain is reached when A^ > A |10|. 

The black curve shown in Fig. 2e show the LDOS prohle extracted along the Fe chain 
at 5 K. The orange curve represents a tentative ht with the theoretical wavefunction 
giving ,^1 110 nm and .^2 ~ 0.75 nm. We remark that .^25 constituting the short localiza¬ 

tion of the MBS, is about the same as the diameter of the halo observed by AFM in Fig. 2c. 
However, the second localization length ,^1 is only indicative as it exceeds the chain length 
and is presumably affected by the defect states present in this chain. 

To improve on such issues, we systematically investigated defect-free chains hosting MBS 
identihed by both AFM and conductance measurements. Figures 3a and b show the zero- 
bias constant-height AFM maps and the normalized LDOS{x,y) of such defect-free chains 
below Tc. Both data were acquired at slightly larger tip-sample separations (~ -|-50 pm) 
compared to Fig. 2c. For that reason, the AFM image do not resolve the atomic lattice 
along the chain. However, both the dark halo (Fig. 3a) and a ZBP in the conductance map 
Fig. 3b remains present at the chain end testifying the MBS. As previously discussed, a non- 
negligible LDOS signal along the chain is detected in the LDOS map Fig. 3b and decays with 
respect to the chain end. We extracted the period of the oscillations by performing a Fourier 
transformation (FT) of the conductance map (See Supp. Infos. Fig. S7), which gives us the 
Fermi wavevector kp equal to 8.3 nm~^, i.e. a Fermi wavelength of 0.76 nm. The data was 
htted with Eq.(2) by keeping constant the experimentally deduced parameters kp. Fig. 3c 
shows the experimental LDOS{x) prohle extracted along the chain (black dots). As shown 
by the superimposed blue curve, a remarkable agreement with the raw data periodicity is 
obtained allowing the extraction of the two localization lengths and ^ 2 - These localization 
lengths are then found to be 22 nm and 0.75 nm, respectively, and correspond to about 59 
and 2 atomic sites with respect to the chain lattice. Several data sets were analyzed this 
way always showing the same values for these localization (see Fig. S5). In addition, the 
proximity gap associated with ,^2 is found to be A^ ~ 33 meV. 

To account for the effect of the tip on the wavefunction measurement, we modify the 
formula for the probability density by including a broadening effect resulting from a 
tip size effect (see additional text in Supp. Info.). This approximation considers a symmet- 


9 


ric apex having metallic character and giving the modulus squared of the wavefunctions at 
the Fermy energy. The orange curve of Fig. 3b shows the result of such an approximation 
using a tip apex of 0.17 nm and demonstrates that such a broadening effect is sufficient to 
reasonably reconstruct the experimental data. We think that the use of lower measurement 
temperatures and p-wave STM tips [29] in future experiments might help to resolve more 
precisely the MBSs wavefunction. 

In conclusion, our results validate the existence of MBSs at the end of atomic iron chains 
on superconducting lead. Atomic-scale AFM imaging shows that Fe adatoms form mono- 
atomic and straight chains on the Pb(llO) surface, where zero-bias conductance peaks 
emerge at their extremities which we interpret as signature for MBS. These conductance 
peaks do not survive the suppression of superconductivity when increasing the sample tem¬ 
perature, in accordance with the expected behavior of MBSs. Comparison between current 
and force channels further demonstrates that the AFM imaging is sensitive to the MBS 
observed as an additional force contribution. Furthermore, we spatially characterized the 
localization of the MBS wavefunction which is composed of two localization lengths extend¬ 
ing up to about 60 atomic sites. This suggests that the relatively short localization length 
of the MBS in such hybrid Fe/Pb systems arises from weak magnetic coupling of the chain 
atoms and might be tuned by weak external magnetic fields [10, 17]. 

Method 

Sample Preparation Single Pb(llO) crystal were provided by Mateck GmbH. After an 
ex-situ chemical treatment with hydrogen peroxide and acetic acid solution, the sample were 
atomically cleaned under ultrahigh vacuum by several cycles of sputtering and annealing. 
Iron atoms were deposited from a heated e-beam evaporator with a rate of ~ 0.07 ML/min 
on the surface and annealed at ~ 400 K in order to promote the formation of chains. 

Scanning Probe Microscopy. The SPM measurements were realized with a low- 
temperature STM/AFM microscope (Omicron Nanotechnology GmbH) based on a tuning 
fork sensor (/o ~ 25 kHz, k ~ 1800 N/m) and operated at ~ 5 K in ultrahigh vacuum 
(UHV). All STM images were recorded in the constant current mode with the bias voltage 
applied to a bulk tungsten tip. The very-end of the tip apexes were prepared into a 
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clean Cu(lll) surface by gentle indentations. Conductance measurements were done at 
constant-height and zero-bias using the lock-in technique (/ = 570 Hz, A = 200 /iV) with 
the feedback loop opened at / = 100 pA, V = 10 mV. AFM images were performed in the 
constant-height mode at zero-bias with oscillation amplitudes of A = 50 pm. The AFM 
images were obtained in the constant-height mode with a tuning fork sensor oscillating 
with amplitude of 50 pm. The variation of the resonance frequency, A/, resulting from 
site-dependent interaction forces between tip and sample are dynamically tracked using a 
phase-lock loop while the probe travels over the surface. 
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